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SIMPLIFIED COMPOSITE MICROMECHANICS FOR PREDICTING MICROSTRESSES 


Christos C. Chamls* 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

A unified set of composite micromechanics equations Is summarized and 
described. This unified set Is for predicting the ply microstresses when the 
ply stresses are known. The set consists of equations of simple form for pre- 
dicting three-dimensional stresses (six each) In the matrix, fiber, and Interface. 
Several numerical examples are Included to Illustrate use and computational 
effectiveness of the equations In this unified set. Numerical results from 
these examples are discussed with respect to their significance on microcrack 
formation and, therefore, damage Initiation In fiber composites. 


1. INTRODUCTION 

The stresses (strains) (type and magnitude) In the matrix, In the fiber, 
and at the fiber/matrix Interface, determine the structural Integrity and/or 
the damage tolerance of fiber composites. These stresses are normally called 
microstresses In the fiber composites community. Microstresses are present In 
fiber composites from the Instant the composite starts curing. Other micro- 
stresses are Induced by the applied loads and/or moisture pick-up. In general, 
there are three kinds of microstresses: hygral (moisture), thermal, and 

mechanical. These microstresses vary In magnitude, sense (tension and com- 
pression), and type (normal and shear). At each point In the fiber. In the 
matrix, and/or at the Interface there generally are six components of stress 
(strain) - three normal and three shear. 

Determination of the magnitude, sense and type of microstresses requires 
highly Intricate experiments (such as photoelasticity) and/or very detailed 
three-dimensional finite element analysis. In few limited cases, the mag- 
nitude of transverse and shear microstresses can be estimated by simple theo- 
retical expressions. These simple expressions are used to determine the stress 
(strain) magnification between fibers and, therefore, obtain an assessment of 
possible microcracking and Interfaclal disbonds. Recent research at NASA Lewis 
focused on deriving simple theoretical expressions for all six microstresses 
In the fibers, In the matrix, and at the Interface. As a result, a unified set 
of simplified micromechanics equations was developed to predict the fiber com- 
posite microstresses. The objective of this report Is to describe this unified 
set and Illustrate Its significance through appropriate numerical examples. 


♦Senior Research Engineer, Aerospace Structures/Composites. 



2. DEFINITIONS AND CONSTITUENT MATERIALS 


The branch of composite mechanics which provides the formal structure 
whereby the ply microstresses are related to the ply stresses and to the ply 
constituent materials Is called composite micromechanics. Composite micro- 
mechanics for predicting ply microstresses Is defined concisely In figure 1 
which also Identifies the required Input and the resulting output. The Input 
consists of: (1) constituent materials and their respective properties, 

(2) ply geometric configuration, (3) fabrication process variables, and 
(4) ply stresses. The output consists of: (1) ply microstress Influence 

coefficients, (2) fiber stresses, (3) Interfaclal stresses, (4) Interfiber 
matrix stresses, (5) Intermatrix matrix stresses, (6) fracture checks for local 
fiber, Interfaclal bond, and local matrix, and (7) ply microcracking assessment. 

The formal structure of composite micromechanics (concepts, math-models, 
and equations) for predicting ply microstresses Is derived based on assump- 
tions which are consistent with (1) the physical situation and (2) the funda- 
mental principles of solid mechanics. These assumptions are: (1) the ply 

resists primary loads by In-plane action as depicted schematically In figure 2, 
(2) the ply resists secondary loads by through-the-thlckness (Interlaminar) 
action, and (3) the ply and Its constituents behave linearly elastic, or nearly 
so, to fracture as depicted schematically In figure 3. Though the principles 
of solid mechanics can be used with various levels of sophistication, the ele- 
mentary mechanlcs-of -material s theory was used In deriving the equations sum- 
marized herein because: (1) It leads to explicit equations of simple form and 

(2) It Is consistent with the simplified composite micromechanics equations 
used to predict ply hygrothermomechanlcal properties (ref. 1) and ply 
strengths (ref. 2). 

The ply microstresses are determined through composite stress progressive 
decomposition as depicted schematically In figure 4. The ply stress can be 
determined from composite stresses as described In references 4 to 7. The 
composite stresses (o c ) are defined In the composite structural coordinate 
reference axes (x,y,z). The ply stresses ( o®.) are defined In the ply 
material coordinate reference axes (refs. 1, 2, 3). The ply microstresses are 
defined In a representative Intraply region In the ply material coordinate 
reference axes In two subregions. Interfiber (B) and intermatrix (A), as shown 
at the lower left In figure 4 and. In more detail In figure 5. The Interfiber 
microstresses In the matrix In subregion (B) are denoted by (fig. 5) and 

those In the fiber by (fig. 5). The Intermatrix micro-stresses In the 

matrix In subregion (A) are denoted by (fig. 5). The Interfaclal 

microstresses are assumed to be: (1) "zero" In subregion (A), (2) the same as 

/ D \ 

those In the Interfiber matrix In subregion (B) but denoted by , and (3) 

combinations of subregions (A) and (B) Identified as subregion (C) In figure 6 
(C) 

and denoted by a ^ N 

Typical properties of select fibers, and matrices, and composites are 
respectively summarized In tables 1 and 2. Properties from these tables are 
used to determine ply microstresses as described In the Illustrative examples. 
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3. PLY MICROSTRESSES DUE TO PLY LONGITUDINAL STRESS (o all ) 

The equations to predict stresses In the constituents due to a uniaxial 
longitudinal ply stress are summarized In figure 7. The schematics In the 
figure define the notation In the equations and their respective subregions. 

The longitudinal matrix microstresses are assumed to be the same In both sub- 
regions (A) and (B). The longitudinal microstresses depend on relative con- 
stituent and ply moduli while the transverse microstresses also depend on 
constituent and ply Poisson's ratios and expllcltely on the fiber volume ratio. 

In order to use the ply microstress equations In figure 7, we need two 
fiber properties (Efn and »fi 2 )* two matrix properties (E m and v m ) , 
two ply properties (Em, vjtl 2 )» the fiber volume ratio ( kf ) and the ply 
longitudinal stress (qn). The fiber and matrix properties are obtained from 
tables 1 and 2. The ply properties are calculated from the equations listed 
at the bottom of figure 7. The following examples Illustrate use of the equa- 
tions. 

Example 3.1.: Determine the constituent stresses (ply microstresses) In 
an AS/IMHS ply with 0.6 fiber volume ratio, assuming: (a) 200 ksl tensile ply 
stress, (b) 180 ksl compressive ply stress and (c) <^23 = 0. These ply 
stresses correspond to ply uniaxial fracture stresses. 

Solution: The required computations are expedited by proceeding as fol- 
lows: ( 1 ) select respective constituent material properties from tables 1 and 2 

(2) calculate ply properties from the equations In figure 7; (3) substitute 
these properties In the microstress equations together with the given fiber 
volume ratio and ply stress(es); (4) calculate the microstresses from the 
appropriate equations In figures 7 and 6 (for Interfaclal bond); and (5) sum- 
marize the results. For this example, the constituent material and ply prop- 
erties are: 


Property, units 


Modulus, mpsl 
Poisson's ratio 


Property type 

Fiber 

Matrix Ply 

31 .0 

0.5 18.8 

0.20 

0.35 0.26 


The corresponding ply microstresses (rounded to one decimal place) are: 


Subregion Ply microstresses (ksl) for ply stress 


mlc 

rostress 

o a -|-| = 200 ksl 

° 8,11 = _18 ° 

°mll 

(A or B) 

5.3 

-4.8 



330.0 

-297.0 

(A) 
°m 22 * 

(A) 

°m33 

0.5 

-0.4 

JB> 

°m 22 ’ 

JB) 

°m33 

- 0.1 

0.1 

JB) 
°f 22 ’ 

JB) 
°f 33 

- 0.1 

0.1 
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Ply microstresses (ksl) for ply stress 
°#,n = 200 ksi = ~ibo ksi 

-0.1 o.i 


0.6 -0.4 

0 0 

Examining the magnitudes of the microstresses, the following are observed: 
( 1 ) the longitudinal ply stresses are primarily resisted by the fibers as 
expected; (2) The stresses In the matrix In subregions (A) and (B) are negli- 
gible compared to matrix strength (15 ksl tensile, 35 ksl compressive In 
table 2); (3) ply longitudinal stress Induces microstresses In subregions (A) 
and (B) of opposite sign; and (4) the Interfaclal bond stress In subregion (C) 
may be of sufficient tensile magnitude to cause Interfaclal disbonding. (The 
Interfaclal tensile bonding strength Is typically estimated to be about 1.0 ksl). 
It may be concluded, therefore, that longitudinal ply stresses will generally 
Induce microfractures associated with (In the order of most probable occurr- 
ence): (1) fiber fractures, (2) Interfaclal disbonds, and (3) matrix fracture. 


Subregion 

microstress 

(B) (B) 

*122’ °133 

, (c) 

J 1N 

,(C> 

IS 


4. PLY MICROSTRESSES DUE TO PLY TRANSVERSE STRESS (<tj.22) 

The equations to predict stresses In the constituents due to a transverse 
ply stress are summarized In figure 8 . The schematics In the figure depict the 
notation In the equations and their respective subregions. The longitudinal 
microstresses In the matrix are the same In both subregions (A) and (B). These 
microstresses depend on relative constituent/ply Poisson's ratios and moduli. 
The transverse microstresses depend on relative local/ply moduli (ply moduli 
defined at the bottom of the figure). These microstresses also depend on the 
constituent properties and the fiber volume ratio through the ply moduli and 
Poisson's ratio. 

In order to use the ply microstress equations In figure 8 , we need four 
fiber properties (Ef^, Ef 22 t «f! 2 * and w f23)» two matrix properties 
( Em and w m ), four ply properties (Em, Ejt 22 . * 112 . and »i23). the fiber 
volume ratio (kf), and the ply transverse stress (<j& 22)- The fiber and 
matrix properties are obtained from tables 1 and 2. The ply properties are 
calculated from the equations shown at the bottom of figures 7 and 8 . The ply 
fiber volume ratio and the ply transverse stress are assumed known. The fol- 
lowing example Illustrates use of these equations. 

Example 4.1: Determine the constituent stresses (ply microstresses) In 

an AS/IMHS ply with 0.6 fiber volume ratio assuming: (a) 7 ksl tensile trans 

verse ply stress, (b) 35 ksl compressive transverse ply stress and (\)a ^ 23 = 0. 
These ply transverse stresses are about equal to corresponding uniaxial ply 
fracture stresses. 

Solution: The required computations are expedited by proceeding as fol- 
lows: ( 1 ) select respective constituent material properties from tables 1 

and 2 ; ( 2 ) calculate ply properties from the equations at the bottom of 
figures 7 and 8 ; (3) substitute these properties In the microstress equations 
together with the ply transverse stress; (4) calculate the microstresses from 
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the appropriate equations In figures 8 and 6 (for Interfaclal bond); and 
(5) summarize the results. For this example, the constituent material and ply 
properties are: 


Property, units Property type 



Fiber 

Matrix 

Ply 

Modulus, mpsl (1,1) 

31 .0 

0.5 

18.8 

(2,2) or (3,3) 

2.0 

0.5 

1 .0/1.23 

Poisson's ratio, 




(1.2) or (1,3) 

0.20 

0.35 

0.26 

(2,3) 

0.25 

0.35 

0.42 


The corresponding ply microstresses (rounded to one decimal place) are: 


Subregion 

microstress 

(A or B) 
°ml 1 


Ply microstresses ( ks 1 ) for ply stress 
a %22 = 7 a 9.22 ~ ~ 35 ksl 


2.4 -12.0 


“fll 

(A) 
s m22 
3 (B) 
’m2 2 

,< B > 

3 f 22 

, (A) 

J m33 

(B) 
5 m33 

(B) 
J f 33 
(B) 
J 122 
,(B) 
*133 

, (C) 

*1N 

,< C > 

*1S 


-1.6 

8.0 

3.5 

-17.5 

8.4 

-42.0 

8.4 

-42.0 

-0.2 

1.5 

0.1 

-0.3 

0.1 

-0.3 

8.4 

42.0 

0.1 

-0.3 

8.3 

-41.6 

6.0 

-30.1 


*(^ 22 ^ 122 ^ as P red ^ cted ^ rom the respective equations In figure 8. 
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In examining the magnitudes of the above ply microstresses, the following 
are observed: ( 1 ) the ply transverse stress Is resisted mostly by the con- 

stituents In subregion (B); (2) the longitudinal stress In the matrix has the 
same sign as the ply transverse stress; (3) the longitudinal stress In the 
fiber has the opposite sign of the ply transverse stress and Is of relatively 
negligible magnitude compared to corresponding strength; (4) the through-the- 
thlckness microstresses In subregions (A) and (B) are negligible In both 
constituents; ( 5 ) the normal Interfaclal bond stress Is of relatively high 
magnitude In subregion (B); and ( 6 ) both Interfaclal stresses are of relatively 
high magnitude In subregion (C). It may be concluded, therefore, that the 
transverse ply stress will generally Induce microcracks associated with (In the 
order of most probable occurrence): ( 1 ) Interfaclal bond and ( 2 ) Intermatrix 

matrix fracture. 


5. PLY MICROSTRESSES DUE TO PLY FLAT-WISE STRESS ( 0 * 33 ) 

The ply flat-wise stress Induces microstresses which are similar to those 
Induced by the ply transverse stress but rotated 90°. The similarity Is a 
consequence of the square array unit cell assumption which has the same Inter- 
fiber and Intermatrix structure, the same constituent properties along the 2 
and 3 material axes direction, and the same ply properties along the 2 and 3 
material axes directions. 

The equations to predict stresses In the constituents due to a ply 
"flat-wise" stress are the same as In figure 8 , but with the subscripts 22 and 
33 Interchanged. In order to calculate the ply microstresses, first Inter- 
change the subscripts In the equations In figure 8 and second follow the pro- 
cedure described In Example 4.1. The observations and conclusions for this 
example, relative to microstress magnitudes and probable microfractures, are 
applicable for the ply flat-wise stress as well but, again, rotated 90°. 


6 . PLY MICROSTRESSES DUE TO PLY INTRAPLY SHEAR STRESS (oj,i 2 or <j al3 ) 

The equations to predict stresses In the constituents due to an Intraply 
(Intralaminar) shear stress are summarized In figure 9. The schematics In the 
figure depict the notation used In the equations and their respective sub- 
regions. The microstresses depend on relative constituent/ply shear moduli. 
The dependence of the microstresses on the fiber volume ratio Is through the 
ply shear moduli G-j 2 and as shown In the figure. The equations 

for the shear microstresses In the 13 plane due to ply Intraply shear stress 
<*^13 are the same as those In figure 9 with the subscripts 12 replaced by 
13. 


In order to use the ply microstress equations In figure 9, we need one 
fiber property (Gfj 2 ). one matrix property (G m ), two different ply sheer 
moduli ( G *( 2 ancl Gj,l 2 ) » the fiber volume ratio (kf), and the Intraply 
shear stress ( 0 ^ 12 ) • The fiber and matrix properties are obtained from 
tables 1 and 2. The two ply sheer moduli are calculated from the equations 
shown In the figure. The ply fiber volume ratio and the Intraply shear stress 
are assumed known. The following example Illustrates use of the equations. 

Example 6.1.: Determine the constituent stresses (ply microstresses) In 

an AS/IMHS ply with 0.6 fiber volume ratio and 10 ks 1 Intraply shear stress. 
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An Intraply stress of 10 ksl Is about equal to the Intralaminar fracture 
stress. 

Solution: The required computations are expedited by proceeding as fol- 
lows: (1) select respective constituent material properties from tables 1 

and 2; (2) calculate ply shear moduli from the respective equations In 
figure 9; (3) substitute these properties In the microstress equations together 
with the Intraply shear stress; (4) calculate the microstresses from the appro- 
priate equations In figures 9 and 6 (for the Interfacial bond); and (5) sum- 
marize the results. For this example, the constituent material and ply 
properties are: 


Property, units 


Property type 
Fiber Matrix Ply 


Shear modulus, mpsl 2.0 


0.19 0.54/0.64 3 


The corresponding ply microstresses (rounded to one decimal place) are: 


Ply microstresses (ksl) 
for Intraply shear stress 
<r t i 2 = 10 ksl 

3.5 

11.9 

11.9 

11.9 

0 

Examining the magnitudes of the above microstresses the following are 
observed: (1) the shear microstress In the matrix In the Interfiber region 

(B) Is substantial and about 90 percent of the matrix shear strength (13 ksl, 
table 2); (2) the Interfaclal shear stress exceeds the ply Intralaminar shear 
strength (about 10 ksl); (3) the shear microstress In the matrix In subregion 
(A) Is relatively low compared to both Intralaminar shear strength and matrix 
shear strength. It may be concluded, therefore, that the Intraply shear stress 
will generally Induce microcracks associated with (In the order of most prob- 
able occurrence) (1) Interfaclal bond or Interfiber matrix fracture and 
(2) Intermatrix matrix fracture. Interfiber matrix fracture Is equally prob- 
able since the magnitude of the shear microstress Is so close to the shear 
strength of the matrix. It Is the equally probable occurrence of either 
Interfaclal bond fracture or Interfiber matrix fracture that explains, In part, 
the sensitivity of Intralaminar or Interlaminar shear strength to: (1) Inter- 

faclal bonding conditions, (2) environmental effects on matrix dominated com- 
posite properties, and (3) matrix effects on composite properties, In general. 


Subregion 

microstress 


« (A) 

°ml2 

(B) 

°ml2 

(B) 

°f 1 2 
(B) 

°112 

All others 


a& 12 /G !ll2' 


7 



As was already mentioned, ply Interlaminar shear stress (0^13) will 
Induce the same microstress magnitudes In the 1-3 plane. The calculation 
procedure, observations, and conclusions are also the same. 


7. PLY MICROSTRESSES DUE TO PLY THROUGH-THE-THICKNESS SHEAR STRESS (o a 23) 

The equations to predict the stresses In the constituents due to a through- 
the-thlckness shear stress are summarized In figure 10. The schematics In the 
figure depict the notation used in the equations and their respective sub- 
regions. The microstresses depend on relative constituent/ply shear moduli. 

The dependence of the microstress on the fiber volume ratio Is through the ply 
shear moduli G23 and Gj^ as shown In the figure. 

In order to use the ply microstress equations In figure 10, we need to 
know one fiber property ( Gf 23) » one matrix property (Gm), two ply prop- 
erties (G23 and Gjt23) . the fiber volume ratio (kf), and the through- 
the-thlckness ply shear stress (<Jj,23)- The fiber and matrix properties 
are obtained from tables 1 and 2. The ply properties are calculated from the 
equations shown in the figure. The ply fiber volume ratio and the ply 
through-the-thlckness shear stress are assumed known. The following example 
Illustrates use of the equations. 

Example 7.1.: Determine the stresses In the constituents (ply micro- 

stresses) In an AS/IMHS ply with 0.6 fiber volume ratio assuming a 6 ks 1 
through-the-thlckness ply shear stress. This shear stress Is approximately 
equal to the corresponding ply shear strength which Is estimated to be about 
60 percent of the Intralaminar shear strength (10 ks 1 previous section). 

Solution: The required computations are expedited by proceeding as fol- 

lows: (1) select respective constituent material properties from tables 1 
and 2; (2) calculate the ply shear moduli from the respective equations In 
figure 10; (3) substitute these properties In the microstress equations 
together with the ply through-the-thlckness shear stress; (4) calculate the 
microstresses from the appropriate equations In figures 10 and 6 (for Inter- 
facial bond); and (5) summarize the results. For this example, the con- 
stituent material and ply properties are: 

Property, units Property type 

Fiber Matrix Ply 

Shear modulus, mpsi 1.0 0.19 0.51/0.37 3 


aG 23 /G H23- 
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The corresponding ply microstresses (rounded to one decimal place) are: 

Ply microstresses 
(ksl) for through- 
the-thlckness ply 
shear stress 
0*23 - & ksl 

3.1 

8.3 

8.3 

8.3 

6.0 

0 

Examining the magnitudes of the above microstresses the following are 
observed: (1) the shear microstresses In the matrix In the Interfiber sub- 

region B are substantial compared to the corresponding ply shear strength 
(6 ksl); (2) the Interfaclal bond shear strength In the Interfiber subregion B 
Is also substantial; (3) the Interfaclal bond normal stress In subregion (C) 

Is equal to the through-the-thlckness shear strength; and (4) the matrix 
microstress In the Intermatrix subregion (A) Is about 50 percent of the cor- 
responding ply shear strength. It may be concluded, therefore, that the 
through-the-thlckness ply shear stress will generally Induce microcracks asso- 
ciated with (In the order of most probable occurrence) (1) Interfaclal bond or 
Interfiber matrix fracture and (2) Intermatrix matrix fracture. If It Is 
assumed that the Interfaclal bond Is the weakest link, then the microcracks 
will develop at this location first. Subsequently the microcracks will pro- 
gress to the Intermatrix subregion and produce transply crack or ply fracture. 

It Is Interesting to note that the Interfaclal bond microstress magnitudes 
for this example are almost Identical to those for the ply transverse tensile 
stress In Example 4.1. One conjecture may be that both of these ply stresses 
Induce ply fracture which Initiates at the Interface and, therefore, the cor- 
responding ply strengths are Interfaclal bond-strength controlled. This Is 
consistent with the general consensus view In the composites community. 


8. PLY MICROSTRESSES DUE TO TEMPERATURE CHANGE (AT*) 

The equations to predict the stresses In the constituents due to a tem- 
perature change In the ply are summarized In figure 11. The schematics In the 
figure depict the notation used In the equations and their respective sub- 
regions. The ply microstresses due to temperature change depend on the rela- 
tive constituent materlals/ply thermal expansion coefficients, the constituent 


Subregion 

microstress 


0 (A) 

*m23 

0 (B) 

*m23 

0 (B) 

*f 23 

0 (B) 

*123 

0 (C) 

*1N 

All others 
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materials moduli, and on the fiber volume ratio. These microstresses also 
depend on the fiber volume ratio and the ply longitudinal modulus through the 
ply thermal expansion coefficients as shown In the figure. 

In order to use the ply microstress equations In figure 11, we need to 
know three fiber properties (afy|, <*f22» and E fll)* three matrix prop- 
erties (am, E m and v m ), three ply properties (asm. <*t 22 . and Em), 
the fiber volume ratio (kf), and the temperature change In the ply (ATj,). 

The fiber and matrix properties are obtained from tables 1 and 2. The ply 
properties are calculated from the equations shown In the figure. The ply 
fiber volume ratio and the ply temperature change are assumed known. The fol- 
lowing example Illustrates use of the equations. 

Example 8.1: Determine the stresses In the constituents (ply micro- 

stresses) In an AS/IMHS ply with 0.6 fiber volume ratio and only subjected to 
300 °F ply temperature change. This temperature change corresponds to the 
temperature difference between cure temperature and room temperature. It was 
selected to determine the magnitude of the residual microstresses. 

Solution: The required computations are expedited by proceeding as fol- 

lows: (1) select respective constituent material properties from tables 1 
and 2; (2) calculate the ply longitudinal modulus and the ply thermal expansion 
coefficients from the respective equations In figure 11; (3) substitute appro- 
priate properties In the microstress equations together with fiber volume ratio 
and the temperature change; (4) calculate the microstresses from the appro- 
priate equations In figures 11 and 6 (for Interfaclal bond); and (5) summarize 
the results. For this example the constituent material and ply properties are: 


Property, unit 

Property type 


Fiber 

Matrix 

Ply 

Modulus , mpsl (1,1) 

31 

0.5 

18.8 

Poisson's ratio (1 ,2) 

0.20 

0.35 

0.26 

Thermal expansion 
coefficients , 
v1n/1n/°F 
(1,1) 

-0.55 

36 

-0.16 

(2,2) 

5.6 

36 

15.3 


The corresponding ply microstresses (rounded to one decimal place) are: 


Subregion Ply microstresses 

microstress ( ks 1 ) for a ply 

temperature change 
ATfc = -300 °F 


°mll 


5.4 


3 fll 

, (A) 

m2 2 

3 (B) 

m2 2 


-3.6 


3.1 


-0.9 
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Subregion 

Ply microstresses 

microstress 

(ksl) for a ply 
temperature change 


ATi = -300 °F 

a , , 

5.4 

ml 1 


"fll 

-3.6 

a (A) 

a m22 

3.1 

(B) 

c m 22 

-0.9 

JB> 

°f 22 

-0.9 

(A) 

3.1 

°m33 

(B) 


^33 

-0.9 

(B) 


°f 33 

-0.9 

(B) 

°122 

-0.9 

(B) 


°133 

-0.9 

(C) 

3.1 

°1N 

J C) 

IS 

0 


Examining the magnitudes of the residual microstresses due to cure tem- 
perature (ATj = -300 °F) the following are observed: (1) the longitudinal 

microstresses are relatively low compared to corresponding constituent material 
strengths; (2) the microstresses In the matrix In both subregions (A) and (B) 
are also low; and (3) the normal microstress In the Interfaclal bond (sub- 
region C) may be substantial compared to a typical tensile bond strength of 
about 1 ksl. It may be concluded, therefore, that the residual microstress 
could Induce microcracks at the Interface In subregion (C). These microcracks 
can subsequently extend Into the Intermatrix subregion (A). 


9. PLY MICROSTRESSES DUE TO MOISTURE CHANGE (M t ) 

The equations to predict the stresses In the constituents due to a mois- 
ture change In the ply are summarized In figure 12. The schematics In the 
figure depict the notation used In the equations and their respective sub- 
regions. The microstresses due to moisture change depend on the relative 
constituent materials ply/moisture expansion coefficients, the corresponding 
moduli, and on the fiber volume ratio. These microstresses also depend on the 
fiber volume ratio and the ply moduli through the ply moisture expansion 
coefficients as shown In the figure. 

In order to use the ply microstress equations In figure 12, we need to 
know two fiber properties (Efn and Ef 22 ) * * wo matrix properties (B m 
and E m ) , four ply properties (Bill. 8122 . E #,11 and E t 22 ). the fiber 
volume ratio (kf), and the moisture change In the ply (Mj.). The fiber and 
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matrix properties are obtained from tables 1 and 2. The ply properties are 
calculated from the equations shown In the figure. The ply fiber volume ratio 
and the ply moisture change are assumed known. The following example Illus- 
trates use of the equations. 

Example 9.1: Determine the stresses In the constituents (ply micro- 

stresses stresses) In an AS/1MHS ply with 0.6 fiber volume ratio assuming only 
1.8 wt % ply moisture content. This moisture change corresponds approxi- 
mately to composite moisture saturation under normal exposure conditions. It 
was selected to determine the magnitude of the ply microstresses due to these 
normal moisture exposure conditions. 

Solution: The required computations are expedited by proceeding as fol- 

lows: (1) select respective constituent material properties from tables 1 and 2; 
(2) calculate the ply moduli and the ply moisture expansion coefficients from the 
respective equations In figure 12; (3) substitute appropriate properties In the 
microstress equations together with fiber volume ratio and the moisture change; 

(4) calculate the microstresses from the appropriate equations In figure 12 and 
for Interfaclal bond In figure 6; and (5) summarize the results. For this example 
the constituent material and ply properties are: 


Property, units 

Property type 



Fiber Matrix 

Ply 

Modulus , mpsl (1,1) 

31 0.5 

18.6 

( 2 , 2 ) 

2.0 0.5 

1.2 

Moisture expansion 



coefficients , 



1n/1n/% M 



(1.1) 

0 0.33 

0.0035 

(2,2) 

0 0.33 

0.0806 


The corresponding ply microstresses (rounded to one decimal place) are: 


Subregion Ply microstresses 

microstress ( ks 1 ) for a ply 

moisture change 
M a = 1.8 percent 


tf mll 

-0.3 

°fll 

0.2 

(A) 

%22 

-0.2 

JB> 

%?2 

0.1 

JB) 
°f 22 

0.1 

<A) 

m3 3 

-0.2 

JB> 

m3 3 

0.1 
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Subregion 

Ply microstresses 

microstress 

(ksl) for a ply 
moisture change 
= 1.8 percent 

JB> 

tf f 33 

0.1 

(B) 

°122 

0.1 

JB> 

°133 

0.1 

(C) 

°1N 

-0.1 

a (C) 

°1S 

0 


Examining the magnitudes of the moisture (hygral) microstresses due to a 
moisture change of 1.8 wt %, the following are observed: (1) the longitudinal 

microstresses are relatively negligible compared to corresponding constituent 
material strengths; (2) the microstresses In the matrix In both subregions 

(A) and (B) are also negligible; and (3) the normal microstress In the Inter- 
face bond (subregion (C)) Is negligible compared to a typical tensile bond 
strength of about 1 ksl. It may be concluded, therefore, that the hygral 
microstress would not Induce microcracks at the Interface In either subregion 

(B) or (C). 


10. DISCUSSION 

The several numerical examples presented and discussed Illustrate the 
usefulness and advantage of having a unified set of micromechanics equations 
summarized In figures 6 to 12 for predicting ply microstresses. The examples 
also Illustrate how the various microstresses are Interrelated and dependent 
on constituent and ply properties. In addition they provide detailed and 
quantitative Information of the ply microstresses In the various subregions and 
also possible sites of microcrack formations. Furthermore, the various micro- 
stress equations can be selectively used to conduct parametric studies as well 
as sensitivity analyses to assess the Influence of constituent materials, 
fiber volume ratio, and mechanical and environmental load conditions on ply 
microstresses and possible microcrack formations. 

The two tables summarizing constituent material properties Illustrate the 
amount of data needed for effective use of a unified set of micromechanics 
equations. The data In these tables were compiled from many sources and many 
values are estimates which were Inferred from predicted results and curve fits. 
The data are Included for three main reasons: (1) to Illustrate that the use 

of these micromechanics equations require numerous property data; (2) to bring 
attention to the fact that many of these properties have not been measured and, 
hopefully, stimulate enough Interest to develop experimental methods to meas- 
ure them; and (3) to provide Indicative ranges of properties of both fibers and 
matrices. It cannot be overemphasized that the data should be considered 
dynamic In the sense that they should be continuously modified If better values 
are known or become available. 
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Lastly, the unified set of micromechanics equations described herein for 
microstresses can be used to determine the effects of voids, moisture, and 
temperature on the magnitudes of the microstresses. This Is accomplished by 
first predicting the effects of these variables on matrix and ply properties 
using the procedures described In references 1 and 2. These modified prop- 
erties are then used In the ply microstress equations to predict the resulting 
effect. The equations can also represent mlcrostress/ply stress Influence 
coefficients simply by dividing the desired microstress with Its respective ply 
stress . 


CONCLUSIONS 

A unified set of micromechanics equations of simple form Is summarized and 
described for predicting ply microstresses. The unified set of equations con- 
sists of simplified expressions to predict microstresses In the matrix, give 
the Interface, and In the fiber for hygral, thermal, and mechanical loading 
conditions. Several examples are presented to Illustrate usefulness of these 
equations and attendant significance to the structural Integrity and damage 
tolerance of fiber composites. The availability of this unified set of equa- 
tions makes It possible to explain the origin of microcracking and provide a 
basis for determining the damage Initiation at the fiber composite microscale. 
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APPENDIX - SYMBOLS 


E modulus of elasticity 

6 shear modulus 

k volume ratio 

M moisture, wt % 

T temperature 

x,y,z structural reference axes 

1,2,3 ply material axes 

a thermal expansion coefficient 

B moisture expansion coefficient 

£,e fracture strain, strain 

v Poisson's ratio 

0 stress 
Subscripts: 

C compression property 

f fiber property or stress 

G glass-transition 

1 Interface 

8. ply property or stress 

m matrix property or stress 

N normal to Interface 

S shear, shear tangential to Interface 

T tension 

v void 

W wet 

1,2,3 direction corresponding to ply material axes 
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Superscripts: 

A microstress 
B microstress 
C microstress 


In Intermatrix subregion (A) 

In Interfiber subregion (B) 

In the Interface In Intermatrix region (C) 
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TABLE 2. - MATRIX PROPERTIES 


Property 




Matrix 



Name 

Symbol 

Units 

LM 

IMLS 

IMHS 

HM 

Poly- 
i mi de 

PMR 

Density 

Pm 

lb/in 
10° psi 
10 6 psi 

0.042 

0.046 

0.044 

0.045 

0.044 

0.044 

Modulus 
Shear modulus 

Em 

^m 

0.32 

0.50 

0.50 

0.75 

0.50 

0.47 

Poisson's ratio 



0.43 

0.41 

0.35 

0.35 

0.35 

0.36 

Heat capacity 

c m 

Btu/lb/°F 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

Heat conductivity 

Km 

Btu/hr/ft 2 /°F/in 

1.25 

1.25 

1.25 

1.25 

1.25 

1.25 

Thermal exp. coef. 

a m 

10-6 in/in/°F 

57 

57 

36 

40 

20 

28 

Diffusivity 


10~ 10 in ? /sec 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

Moisture exp. coef. 

e m 

i n / i n / % M 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

Tensile strength 

Sfflt 

ksi 

8 

7 

15 

20 

15 

8 

Compression strength 

Smc 

ksi 

15 

21 

35 

50 

30 

16 

Shear strength 

s ms 

ksi 

8 

70 

13 

15 

13 

8 

Tensile fracture 

£mt 

1n/in(%) 

8.1 

1.4 

2.0 

2.0 

2.0 

2.0 

strain 








Compr. fracture 

£mc 

in/in (%) 

15 

4.2 

5.0 

5.0 

4.0 

3.5 

strain 








Shear fracture 

£ms 

in/ in(%) 

10 

3.20 

3.5 

4.0 

3.5 

5.0 

strain 








Air heat conductivity 

K v 

Btu/hr/ft^/°F/in 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

Glass trans. temp. 

t gd 

°F 

350 

420 

420 

420 

700 

700 

(dry) 









Notes: LM - Low modulus; IMLS - Intermediate modulus low strength; IMHS = Intermediate modulus high 

strength; HM - High modulus. 


Thermal, hyqral , compression and shear properties are estimates only; G m = E m / 2 ( 1 + \> m ). 
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CONCEPTS, MATH-MODELS AND EQUATIONS USED TO PREDICT FIBER COMPOSITE PLY MICROSTRESSES. 
AND MICROCRACKING KNOWING CONSTITUENT MATERIALS/PROPERTIES. GEOMETRIC CONFIGURATION. 


FABRICATION VARIABLES AND MECHANICAL. HYGRAL AND THERMAL PLY STRESSES 



Figure 1. - Simplified composite micromechanics for predicting microstresses-definition. 




Figure 2. - Typical fiber composite geometry. 













Figure 3. - Typical stress-strain behavior of unidirectional fiber 
composites. 





Figure 4. - Ply microstresses through composite stress progressive decomposition. 


CO-85- W31 
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REPRESENTATIVE INTRAPLY REGION 



INTERMATRIX 
"SUB-REGION" (A) / (A) \ 
MATRIX STRESSES (o m j 


INTERFIBER 

"SUB-REGION" (B) ( (B) 
MATRIX STRESSES (a m 



INTRAFIBER 
"SUB-REGION" (B) (B) 
FIBER STRESSES (a f 


C0-85-l?nn 

Figure 5. - Sub-region ply microstresses definition/notation. 


INTERMATRIX (A) 

(A) (A) 

°i22 * °i33 ~ 0 


INTERFIBER (B) 


(B) 

(B) 


°i 22 * 

°m22 


(B) 

(B) 


°i33 " 

0 m33 


(B) 

(B) 


°i!2 = 

°ml2 


INTERMATRIX (C) 

(C) 1 

17 (A) 

(B) \ 

°iN "( 

[A°m22 

+ a m22> 


( (A) 
+ I°m33 

* 


CNJ 

■3* 

4 - 

I 

<ci r. 

( (A) X 

(B) \ 

°is‘ L " 1 

l°m22 + 

°m22/ 


/ (A) (B) \\ 1 

+ (°m33 + °m33/ J x 7 



MATRIX (m) 
FIBER <f) 
INTERFACE (i) 

PLY U) 


CO-86- 18692 


Figure 6. - Ply microstresses - stresses at the interface. 
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DUE TO LONGITUDINAL STRESS (o £U ) 


°mll " 

LU 

e 

LU 

°eu 

°fll - 

(E fll /E ^ll ) 


(A) 



°m22 = 

tE m 

(B) 

/1-Vk7 \ 

(A) 

°m22‘ 

(vs ) 

°m22 

(B) 

/WS\ 

(A) 

Of 22 “ 

'WS ) 

°m22 

E ill * 

E 

+ 

LU 

^2“ 

E m 


v £\2"~ k i y fI2 + k m v m 

(A) (A) (B) (B) (B) 

a m33 = °m22 s °m33 = a m22 ; °f33 



MICROSTRESS REGIONS 

(B) (B) 


Figure 7. - Ply microstresses - stresses in the constituents. 


DUE TO TRANSVERSE LOAD (0 


°mll = v m 


_ v ll2 E n 


sen 


)°i 22 


t n i 

°fll ' ( v fl2 ' v l\2 22 


°m22 = ( E m/ E 22> J e22 

C = ( E -e22/E22)°i22 

°^2 ’ ( E ^ 22 / E 22)°^22 

(A) 

°m33 * ( v m ‘ v ^23) ( E m l E m) a l22 

(B! .^ 1-^7 j < A > 


°m33 ' 


./W?L\ 

/ 

E 22 = 0 " VS) E m + 


°m33 


(B) 

°f33 


(A) 

%33 


VS E n 






: E i22 “ 


1 'M 1 ' E f2 2 


V^23 ‘ k f v f23 + k m (2v m - V t i2 i t22l i ni) 

Figure 8. - Ply microstresses - stresses in the constituents. 


CO-86- 18690 


3 


DUE TO INTRAPLY SHEAR STRESS (o /12 > 


°ml2 = (G m ,G 12’ °/12 
(B ) 

°ml2 = (g ^12 ,g 12 ) 

(B) . 

°fl2 = * G il2' G 12 °£12 



(A) (A) (B) (B) (B) (B) 

°ml3 = °m!2 = °ml3 = °m!2 : °f!3 * °fl2 


r PLY W) 

PLY MATERIAL AXES / FjBER (f) 



I A | B |A| 



MICROSTRESS REGIONS 


Figure 9. - Ply microstresses - stresses in the constituents. 


DUE TO THROUGH-THE-THICKNESS SHEAR 

(A) 

°m23 = (G m /G i23 ) °i23 

(B) 

0 m23 = (g 23 /g ^23 ) °^23 

(B) (B) 

°f23 * °m23 




PLY MATERIAL AXES 



MICROSTRESS REGIONS 


Figure 10. - Ply microstresses - stresses in the constituents. 
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DUE TO TEMPERATURE (ATg) 

°mll " ( a ill " a m) AT i E m 
°fll * ( a ill " a fll) AT i E fll 
a m22 * ( a i22 ' a m) AT i E m 

<bi in, rr\ ,A) ,B) 

°m22* 'V k f) a m22 * 0 f22 

1 -t£ \ (A) 




(B) (Bl 
°f22 * °m22 


HU * 



[ k f a f 1 1 E fll + k m a m E m]/ E 211 
a i22 * yi<f a f22 + ( E -V^ )( ] + k f u m £ fll / E ill) “it 
£ ill * k f E fll + k m E m 


rPLY (i) 


(A) (A) (Bl (Bl (B) (B) 

°m33 = °m22 : °m33 * °m22 * °f33 * °f22 

Figure 11. - Ply microstresses - stresses in the constituents. 



DUE TO MOISTURE (Mj) 

°mll '(Pill ' Pm) M i E m 

°fll “ Pill M i E fll 

°L A 22 ' ih22 ' ^ E m 

(B) (B) / W 5 \ (A) 

°f22 * °m22 ’ j °m22 

Pill 'l 1 ' k f)Pm E m/ E ill 

VS ' (l ' VS ' ,E " 




V^T E i22 + (l "V k f) E m 

Em 


1 -y[T { (1 - E m /£f 22 ) 


0 (A) . 0 (A) . g (B I . 0 (BI . g(B) = g(BI 

m33 m22 ' m33 m22 ’ f33 f22 



Figure 12. - Ply microstresses - stresses in constituents. 
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